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[ABSTRACT] 

Variational Monte Carlo calculations for \H (ground and excited states) and \He are per- 
formed to decipher information on A-nuclear interactions. Appropriate operatorial nuclear and 
A-nuclear correlations have been incorporated to minimize the expectation values of the energies. 
We use the Argonne V\% two-body AA along with the Urbana IX three-body NNN interactions. 
. The study demonstrates that a large part of the splitting energy in \H (0 + — 1 + ) is due to the 

three-body AAA forces. hypernucleus is analyzed using the s-shell results. A binding to nu- 
clear matter is calculated within the variational framework using the Fermi-Hypernetted-Chain 
technique. There is a need to correctly incorporate the three-body ANN correlations for A binding 
£^ . to nuclear matter. 
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§1. Introduction 

The study of the response of a many-body system to a hyperon gives insight into the struc- 
ture of baryon-baryon interactions. The binding energy data of light s-shell hypernuclei provide 
a unique opportunity to know more about the A-nuclear interactions, particularly on their spin- 
dependence. In the past, basically two approaches have been followed. The first one involves 
Brueckner-Hartree calculations using Nijmegen YN potential with and without higher order cor- 
rection to single-particle energies [1,2]. This method uses the large AN^T,N coupling which gives 
considerably lower binding energy for \He. Any attempt to correct this leads to poor agreement 
with the scattering data. The second method is primarily based on reliable variational techniques, 
mostly using simplified AA interactions [3,4]. We follow this approach but use realistic Argonne 
Ui8 AA interaction [5] along with Urbana IX three-body AAA interaction [6, 7]. The phenomenolog- 
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ical approach we follow is consistent with meson-theoretic models as well as available low energy 
scattering data. The AiV^ETV coupling is effectively taken care of by inclusion of the phenomeno- 
logical ANN potential. The ANN potential consists of both the dispersive and two-pion-exchange 
(TPE) kind as employed in previous studies[4]. The hypernuclei considered in this work are \H, 
\H* ( * on \H refers to the excited 1 + state), and \He. 

The interaction parameters which we find based on our s-shell results are later used to make 
estimates of the binding energy of \ 7 0. We also study the A binding to nuclear matter by us- 
ing the Fermi- Hypernetted- Chain (FHNC) technique[4, 8]. This study gives an indication of the 
implications of our s-shell results on heavier hypernuclear systems. 

In section 2 we describe the Hamiltonian used in this work. Section 3 gives the wave function 
and approach. In section 4 we discuss the results and finally in section 5 we give the conclusion 
and comments. 



§2. The Hamiltonian 

The complete hypernuclear Hamiltonian consists of the nuclear Hamiltonian H^~ x and the 
lambda Hamiltonian Ha- The nuclear Hamiltonian -ff^ -1 is given by 

h^- 1 = - E £- v ? + E v « + E v** > 

i—1 i<j i<j<k 

where Vij and Vijk are the two-nucleon NN and three-nucleon NNN potentials, respectively, and 
mj is the mass of the nucleon. 

The two-body NN interaction employed here is the Argonne vis interaction [5]. The first 
fourteen operator components of this model are charge- independent and are an updated version 
of the Argonne w 14 potential[9]. Three additional charge-dependent and one charge-asymmetric 
operators are added along with a complete electromagnetic interaction, containing the Coulomb, 
Darwin-Foldy, vacuum polarization, and magnetic moment terms with finite-size effects. The 
potential has been fit directly to the Nijmegen pp and np scattering data base[10,ll], low-energy 
nn scattering parameters, and deuteron binding energy. For the three-body NNN potential we use 
the Urbana model[6,7] consisting of the two-pion-exchange (TPE) part of Fujita and Miyazawa[12] 
and a repulsive phcnomenological spin-isospin independent term. We have used the Urbana IX 
model[7] of the interaction where the values of the strength parameters are used in conjunction 
with the Argonne u 18 interaction. 

The lambda Hamiltonian Ha is given by 

h 2 A 1 A 1 

#A = -^Vi + E^A + E^A . ( 2 - 2 ) 

i—1 i<j 

where Via and Vij a are the two-body AN and three-body ANN potentials, respectively, and ttia 
is the mass of the A particle. The first terms of Eqs.(2.1) and (2.2) pertain to the total kinetic 
energy of the nucleons and A, respectively. 

The two-body AApotential Van includes a central potential[4] of the same form for the singlet 
and triplet spin states. These have a theoretically reasonable attractive tail due to the TPE in 
accord with Urbana type potentials [13] with spin- and space-exchange terms. 
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V AN = [(V c (r) -V)(l-e+ eP x ) + ±V a <r A ■ a N ]T^(r) , (2.3) 

where V and V„ are the spin-average and spin-dependent strengths, respectively. P x is the Majo- 
rana space-exchange operator, e is the corresponding exchange parameter, V c is the Woods-Saxon 
repulsive core [4], and T v is the one-pion-exchange (OPE) tensor potential shape modified with a 
cutoff. Further details can be found in Ref . [4] . 

In this study we consider potential parameters that are consistent with low energy Ap scat- 
tering data that essentially determine the value of spin-average strength V — 6.15 ± 0.05 MeV 
[4]- 

For hypernuclei with zero-spin core nuclei, such as A He, the major contribution arises from 
the spin-average strength V while the spin component contributes very little. The spin dependence 
V a is assumed to be positive, which is consistent with hypernuclear spins of mass 4 systems. We 
find that for the s-shell hypernuclei (A < 5) the s-state interaction is dominant but the higher 
partial wave interactions, in particular, the p-state, also make a small but significant contribution 
contrary to earlier studies[14]. The importance of the p-state contribution becomes significant due 
to the A-nuclear correlations. 

Studies on hypernuclei have shown that it is necessary to include a three-body ANN interaction 
in the Hamiltonian. We consider phenomenological ANN forces of the dispersive (spin-dependent 
and spin-independent) as well as the TPE kind 15 ), which arise from the suppression of S, A, . . . 
degrees of freedom by the medium, that is, the second nucleon. 

The dispersive kind has a spin dependence that is given by 

Vfi? N (r ijA ) = W T*(r lA ) T^(r jA )[l + U A ■ ft + <?,-)] ■ (2.4) 
The TPE part of the interaction is given by[15] 

W p = ~ C p (n ■ ?j){X iA , X jA }Y(r iA ) Y{r oA ) , (2.5) 
where X^a is the OPE operator given by 

X kA = (a k ■ <7 A ) + S kA {r kA ) TTr(r kA ) (2.6) 

with 

c , x 3(CT fc -r kA ){a A -r kA ) 

S kA {r kA ) = 5 \ a k ■ cr A ) . (2.7) 

r kA 

In Eq.(2.5) {} represents the anticommutator term. Y^r/^) and T^t^a) are the usual Yukawa 
and tensor functions, respectively, with pion mass /i = 0.7 fm _1 . 

The A-nuclear interaction parameters, V, V^, C p , and W arc considered as unknown. These 
are then fitted as a function of the B A values that have been calculated using the s-shell results. 
Taking these values of V, V a , C p , and W we again perform variational calculations that give us 
the final results for \H, A H* , and \He. These results are later used to analyze \ r O and A binding 
to nuclear matter. 

§3. Wave Function and Approach 
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The trial variational wave function we adopt is of the following form: 



**> = 



1+ E ^+E^a+E^ s + E vit 

i<j<k i<j i<j i<j<k 



II fijk 
i<j<k 



%) • (3-1) 



The pair wave function | ty p ) is a symmetrized product of two-body (1 + Uij) and (1 + Ui\) 
correlation operators acting on a Jastrow trial function. This is written as 



' A-l 

i<3 



s n (i + u iA ) 



»=i 



f7jj in Eq.(3.2) is defined as 

^ = E [ II f? 3 k( r ik> r jk) M r v)0. 



(3.2) 



(3.3) 



p=2,6 k^i,j 



with = [1, CTj • CTj, Sy] ® [1, Tj, -Tj]. 5y is the tensor operator, tr and r are the spin and isospin 
operators, respectively. The factor ff- k suppresses spin-isospin correlations between two nucleons 
in the presence of a third one. 

In principle, the Ui\ correlation will consist of a spin and a Majorana space-exchange operator 

[16]: 

U iA = a cr u a (r llv ) a a • a l + a px u px (r iA )P x , (3.4) 

where a a and a px are variational parameters and P x is the space-exchange operator. In our 
calculations we have not included the space-exchange correlations since the calculations become 
complicated and time consuming. In any case the effect of these correlations is expected to be 
small for s-shell hypernuclei. 

The spin dependent correlation u a given in Eq.(3.4) is defined as 



(./• s A ~/ t A ) 



(3.5) 



where is the spin-average correlation function, and / t A are the solutions of the quenched 
A7V potential in singlet and triplet states, respectively, which are given by the following relation: 

h 2 



V z + V s/t (r AN ) + V£ N 



(3.6) 



where V s / t is the quenched AN potential in singlet/triplet state. [i\n is the reduced mass of the 
A-iVpair, while V£ N is an auxiliary potential. 
The Jastrow wave function | tyj) is given by 



A-l 



A-l 



*./>= [n/c A (^)n/ c M] i*). 



(3.7) 



i=l 



i<j 



Here | $) is an antisymmetric product of single particle wave function with the desired (J,T). The 
initial uncorrelated state $ has no coordinate dependence and is real. For example, consider the 
following $ states for \H and \H* expressed in the spin-isospin basis with the appropriate (J,T) 
states, 
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\H (J=0, T=|) = C U,i,0; i-|,0 3 H* A, 5 + 



2 J 2'"' 2' 2' J j = 

|,iO;-|, |, oVh7 = *i A* (3.8) 



and 



3- 



\H* (J= 1, T=i) = C ( 1, 1, 1; ±, i, 1 ) 3 Hj =i A I , (3.9) 



where C represents the Clebsch Gordon coefhcients[17]. 3 H 2 _! and 3 H._ 2 ! are the uncorrelated 

1— 2 3—2 

_1 1 

<f>s for triton, whereas A i 2 and A f are the spin-down and spin-up states, respectively of the A 

2 2 

particle. 

The spin-orbit correlation U^ s is given by 

U^ s = [u ls (n 3 ) + u lsT (n 3 ) n ■ f ] (L • S) tf . (3.10) 

The eight radial functions, f c (fij), u p =2,%(ri ), u; s (r,j), and ui ST (rij) are obtained from ap- 
proximate two-body Euler-Lagrange equations with variational parameters [18]. 

Ujj™ is a three-body correlation induced by the three-nucleon interaction, Vijk- The other 
correlations incorporated in the wave function are a spatial three-body NNN correlation f9- k , along 
with Uijk that consists of a spin-orbit and an isospin three-body correlation. Further details can 
be found in Ref. [19]. 

The three-body ANN correlation Uij\ has the following form: 

U ij A = V ijA (S 1 ,S 2 ) , (3.11) 

where V^a differs from VijA through the cutoff factor c of the usual Yukawa and tensor functions. 
5\ and 62 are variational parameters which multiply the AAW interaction parameters, C p and W , 
respectively. 

No attempt has been made to vary the two-body AW and three-body NNN correlation param- 
eters^] as their effect has been found to be small[20] and only the variational parameters of the 
wave function pertaining to A have been varied to obtain a minimum in the energy. The optimum 
values of these parameters which are used in our final calculations are given in Tables [1] and [2]. 

We calculate energy expectation values using Monte Carlo (MC) integration[21,22]. The 
expectation values are sampled both in configuration space and in the order of operators in the 
wave function by following a Metropolis random walk[23]. The mathematical expressions used to 
evaluate the energy expectation values are given below. 

The energy expectation value for the pure nucleus is given by 



{En } = (K- 1 1 ' ( } 



where ^> N 1 is the wave function of the mass (A-l) nucleus and H N 1 is the nuclear Hamiltonian. 
The energy expectation value for the hypernucleus is given by 

{Eh) ~ (K 1 n) ' (3 ' 13) 



where ^jj is the wave function of the mass 'A' hypcrnucleus and is the hypernuclear Hamil- 
tonian. 

Therefore, binding energy of A to the hypernucleus is given by 

-B A = (E%) . (3.14) 

The nuclear and hypernuclear wave functions, V I / ^ _1 and ^jj are optimized with respect to 
the variational parameters to obtain the minimum in the energies. 

The B\ value for each hypernuclei is calculated from the variational results using Eq.(3.14). 
The Ba value is thus written as a function of the adjustable parameters in the A Hamiltonian 
Ha, and is used to determine the set of parameters which are consistent with the experimental Ba 
values [26,27]. 

§4. Results and Discussion 

Table[3] gives the variational results for the nuclei, namely, 4 He and 3 H calculated using the 
two-body NN Argonne vis interaction [5] and three-body NNN Urbana IX interaction[6,7] with 
relevant correlations. The numbers appearing in parentheses in all the tables in this work indicate 
the statistical error in the last digit. These calculations have been performed on similar lines as 
those by Wiringa et al. [5,19] and the results conform to theirs. These results also check very well 
with the recent calculations of Forest et al. [24], who use a truncated version of the Argonne Vis 
interaction. 

Next we calculate the energy expectation values for the s-shell hypernuclei, namely \H, \H* , 
and \He using the two-body AW Argonne vis and three-body NNN Urbana IX interactions along 
with the two-body A N and three-body A AW interactions with appropriate correlations incorporated 
in the wave function. Variational calculations have been performed for different values of spin- 
average potential strength, V (6.10, 6.15, and 6.20 MeV). The different values of the space-exchange 
parameter e used in this work are 0.24 (V=6.20 MeV), 0.19 (V"=6.15 MeV), and 0.14 (V"=6.10 
MeV) [25]. 

The bulk calculations consist of the energy expectation values for each V as a function of the 
interaction parameters, V a , C p , and W . In Table [4] we illustrate one such set of results for the 
ground state of \H. Our results demonstrate that the Ba values for \He, \H, and \H* show 
similar trends with the spin-average potential strength V and ANN interaction parameters, C p and 
W a . As expected, Ba increases with V. Ba also increases significantly with the increase in Cp, 
while it decreases with W a . As expected, the dependence on V, C p , and W is more pronounced 
for A = 5 than for A = 4 systems. This result is in accord with the earlier calculations where only 
simplified NN interactions have been used[4]. 

An important goal of the present study is to learn about the role of V a , C p , and W through 
the + — 1 + energy splitting in \H and \H* . We place limits on the values of these parameters, 
consistent with the following experimental Ba values: 

Ba{\H) = 2.22 ± 0.04MeV, 
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B A {\H*) = 1.12±0.06MeV, 
B A ( 5 A He) = 3.12 ± 0.02MeV. 



(4.1) 



Values of B A {\H) and Ba are averages for those of \H and \^ e - Limits on the parameters 

V a , C p , and W are determined by the uncertainties in the experimental Ba values. 

For a given value of V we did a x 2 n t for the calculated energy expectation values according 
to the relation 



where B° K is the corresponding value of Ba for V a = C p = W = for each hypernuclear species. 
The coefficients yi-e are varied to give a minimum in the x 2 that is defined as 



Here ./V is the total number of energy calculations for a particular hypernucleus with different 
values of V a , C p , and W Q . Ba is the calculated value of the A separation energy and A£?a is the 
corresponding Monte Carlo statistical error. The values of the coefficients yi_6, as determined 
from this procedure, are displayed in Table [5]. In all the cases considered the x 2 values arc < 1 
which demonstrates the goodness of the fit. There shall be correlated error bars on the coefficients 
2/1-6 which would be reflected in the uncertainties in determining V ai C p , and W a . We find it more 
convenient to consider the uncertainties in the experimental values while placing limits on V a , C p , 
and W - We hope to compensate some of the uncertainties associated with the y's by giving a 
generous allowance to the experimental ABa values as well as by taking into account the Monte 
Carlo statistical errors in the calculation of the energies. 

We use the coefficients yi-e of Table[5] to obtain a fit with respect to the experimental Ba 
values, treating V a , C p , and W as parameters to determine the best fit. We again construct a x 2 
fit using Eq.(4.3) but now "iV" refers to the factor "3" for the three hypernuclear species and Ba 
refers to the experimental Ba values. The x 2 fit is minimized with respect to C p and W a for a 
given value of V a . In Fig. 1 we plot x 2 , C p , and W as a function of V a for V — 6.15 MeV. It is 
seen that both C p and W decrease with increase in V a , the effect being more pronounced for C p . 
Fig. 2 displays the calculated values of Ba as a function of V a for A H, \H* , and \He. Within the 
accuracy of the graphs the Ba values for \H and \He do not show any dependence on V a in the 
range 0.09—0.26 MeV. As one may expect, the Ba values for \H* depend sensitively on V„, and 
thus in turn on the spin dependence of C p and W . 

The x 2 stays very close to zero (which corresponds to almost an exact fit to Ba (exp) values) for 
V„ = 0.176±0.015 MeV, C p = 1.64^.03 MeV and W Q = 0.026=f0.001 MeV. Most of the deviation 
from zero of the x 2 values in Fig. 1 arise from \H* . The dotted horizontal lines of Fig. 2 display 
the limits on the experimental Ba value of \H* , consistent with the experimental error bar of ± 
0.06 MeV. This places limits on the values of V a , thus, in turn on C p and W . However, the actual 
error bars on these parameters would be larger due to Monte Carlo statistical errors. To take this 
into account, we made a number of energy calculations for V a in the range 0.10 to 0.24 MeV. The 
corresponding values of C p and W Q have been taken from the calculations of Fig. 1. We could 



B A (V a , W , C p ) = Vl V a + y 2 W + y 3 C p + y 4 W 2 Q +y 5 C^+ y 6 W a C p + B 



•A ! 



(4.2) 




-I I 



(4.3) 
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obtain acceptable fits to the energies for V a in the range 0.12—0.23 MeV. This gives for V=6.15 
McV: 

V a = 0.176 ± 0.05, C p = 1.64 t 0.15, W Q = 0.026 =F 0.003. (4.4) 
A similar study for V=6.20 MeV gives 

V a = 0.125 ± 0.05, C p = 1.52 t 0.15, W D = 0.025 =F 0.003. (4.5) 

We have been able to obtain good fits only for V' = 6.15 and 6.20 MeV. The value of V = 6.10 
MeV does not reproduce the correct binding energies of s-shell hypcrnuclci. Therefore, we have 
not carried out any error analysis for t^=6.10 MeV. For the sake of completeness wc mention the 
best parameter values for V=6.W MeV: 

V a = 0.193, Cp = 1.84, W = 0.027. (4.6) 

We can note from Table[6(A)] that the AA^ spin potential has a non zero contribution even in a 
closed-shell system such as \He. This arises because of the AAf spin-spin correlations incorporated 
in the wave function. 

Comparing the results for the core nuclei (Table [3]) with the results for the hypcrnuclei 
(Tables[6(A-C)]) wc note, in general, a shrinking of the core nuclei by about 20% in all the hy- 
pernuclei due to the presence of the A particle. This decrease in radii of the core nuclei would 
imply that the lambda wave functions are closer for larger A. This also contributes to the fact that 
the dependence of B\ on V, C p , and W is more pronounced for the mass 5 than for the mass 4 
hypernuclei. The change in the <i-state probability is found to be small in all cases. 

Table[7] gives the breakup of the + — 1 + splitting contributions in A H and \H* arising from 
the AAT spin-dependent strength V a and the three-body ANN interaction Van n for different values 
of spin-average strength V. It can be noted from this table that the energy difference between \H 
and \H* is consistent with the total contribution from V„ and Vann within the error bars of the 
MC calculations. It can be seen that a large part (~ |) of the splitting comes from the three-body 
AAA potential. The two-body contribution arising from V a is around ~ | of the total splitting. 
This is in contrast to the earlier studies [4,28,29] wherein the + — 1 + splitting has been thought 
to have arisen mainly from the spin dependence of the two-body A AT potential. The present study 
clearly demonstrates that Vann plays a significant role in explaining the splitting. This also results 
in a reduced V a as compared to the value of 0.23±0.02 found in Ref. [4], though in our case the 
error bar on V„ is much larger due to reasons discussed earlier. In the present study for V^-=0.23, 
half of the splitting arises because of V a and the remaining half from the three-body ANN forces. 
For the extreme case, in particular, for V a =0.12, the three-body forces contribute nearly | of the 
total splitting. It would be desirable to have an independent fix on V a , for example, from a more 
refined Ap scattering data. This can enlighten us further on the three-body AAA forces. 

The Majorana space-exchange contribution for the various hypernuclei have been found to be 
small but significant in all s-shell hypernuclei. It is in the range 0.1 — 0.7 MeV for e in the range 
0.14 — 0.24 and as expected has a linear dependence with the Majorana exchange parameter e. 

A few variational calculations have been carried out with only the ANN potentials and no ANN 
correlations. Wc find that without the A AW correlations \ He is not bound. We also notice that the 
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contributions from Va N n and V^ N become more repulsive without the correlations and the total 
contribution from the three-body ANN potentials become positive thereby decreasing Ba- We have 
performed a few energy calculations for \He, A H, and \H* but with the three-body AAA part of 
the Hamiltonian completely switched off . We find in this case that \He is overbound by about 
2.34 MeV for V = 6.20 MeV and by about 1.39 MeV for V = 6.15 MeV. In general, the results for 
\H show that it is underbound while A H* is overbound without the three-body ANN interactions 
for all values of V. Both these studies show the importance of the three-body ANN potentials and 
correlations in obtaining a consistent fit to the Ba values for all the s-shell hypernuclei considered 
here. In particular, we notice the importance of the ANN correlations fW N on the effect of TPE 
ANN forces V^ N . These correlations reduce reasonably the repulsive three-body contribution to 
an attractive contribution implying a strong non-linear dependence on C p . 

Implications of s- shell results on \ r O: We now examine the A 7 hypernucleus in relation to the two- 
and three-body A-nuclear potential parameters that we find from our analysis of s-shell hypernuclei 
as described earlier. Usmani, Pieper, and Usmani 16 ) (referred to as UPU) have carried out MC 
calculations for the \ 7 hypernucleus using the vq part of the older Argonne t>i4 potential [9]. 
For the three-nucleon potential they use the same form (Urbana model) as in the present study 
but with different strength parameters (A„=— 0.0333 and £/ o =0.0038). Their trial wave function 
consists of pair and triplet operators acting on a single-particle determinant. In many respects, 
these calculations are similar to our present calculations. The difference lies in the treatment 
of non-central correlations for which they use the Cluster Monte Carlo method with up to four- 
baryon clusters. The central correlations are treated exactly. We carry out this study in the hope 
of analyzing further our estimated s-shell A-nuclear parameters. UPU have given the following 
empirical relation for C p in the range 0-1 MeV and W in the range 0-0.02 MeV: 

B A = 27.3- 8.9 C p + 11.2 C 2 V +870.0 W Q . (4.7) 

This equation relates the B\ of \ 7 with C p and W a . In order to test the consistency of our 
results with Ba{\ 7 0), we assume that relation (4.7) holds for our values of C p and W . UPU have 
done calculations for spin-average strength U=6.16 MeV with space-exchange parameter e= 0.3. 
Our values of U=6.15 MeV with e=0.17 are closest to theirs. We thus need to modify (4.7) for 
our values of V and e. Unfortunately, there is no simple method to scale relation (4.7) for U=6.15 
MeV, as the scaling can be considerably non-linear. However, since the two values of V are very 
close, we assume that this will not affect the results much. The correction for e is simple, since 
in the absence of space-exchange correlation the space-exchange energies are expected to be linear 
with e. Thus, relation (4.7) can be modified as: 

B A = 27.3 + (e - e a )( (1 - P x ) Van) - 8.9 C p + 11.2 C 2 p + 870.0 W a , (4.8) 

where e o =0.3 as taken by UPU, P x is the space-exchange operator and (Van) is the energy expec- 
tation value of the AA potential. 

Using the entries for v (r)(l — e) and v (r) eP x from Table II of Ref. [16] and using our values 
of C p = 1.6407 and W a = 0.0255 for V = 6.15 MeV we obtain 

Ba = 23.3 ± 1.6 MeV . 
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This is considerably larger than the empirical estimate 16 ' ~ 13.0 ± 0.4 MeV. Thus the results of 
the present study are incompatible with those of Ref. [16]. The reason for this incompatibility may 
largely lie in the use of vq part of the vig hamiltonian for ^O. Another reason can be attributed 
to the use of relations (4.7) and (4.8) for large values of C p and W a for which these relations may 
not be adequate. This discrepancy can probably be resolved by carrying out calculations for \ 7 
with Argonne vig Hamiltonian, which, at the moment is an extremely challenging task. 

Implications of s-shell results on A binding to nuclear matter. The presence of a A particle inside 
nuclear matter can reveal information on the A-nuclear interactions. The well depth D is identified 
with the separation energy for a A in nuclear matter. It is an important parameter which can 
help to distinguish between different AN potentials and also throw light on the AAA interaction. 
A binding to nuclear matter can put further constraints on the potential parameters, namely V, 
V a , W , and C p . With this aim in mind we have performed calculations for D using the Fermi- 
Hypernetted- Chain (FHNC) techniquc[4] to calculate the energy expectation values. 

We have calculated the well depth D variationally using the same underlying principle as 
for our s-shell hypernuclei. Our discussion on D is based on the results given in Table[8]. The 
empirical value of D is now fairly well established at 29 ± 1 MeV [30,31] at the normal nuclear 
matter density of p Q w 0.16 fm~ 3 . These results clearly indicate that A is underbound at the normal 
nuclear matter density, p « 0.16 fm~ 3 . This indeed is a disturbing feature. Bodmer and Usmani[4] 
have found that it is possible to obtain a consistent phenomenology with hypernuclear interactions 
which include the s-shell and the medium and heavy hypernuclei as well as A binding to nuclear 
matter. Our results indicate that with the present available techniques of treating the nuclear 
matter this is not possible. The resolution of this paradox perhaps lies in the proper handling of 
the three-body correlations, particularly the AAA correlations for nuclear matter. It may be noted 
that the contribution from the TPE AAA forces Vfijj N for nuclear matter is always positivc[4]. 
On the other hand, V^ N for s-shell hypernuclei and \ 7 C> is always negative and substantial. The 
three-body AAA correlations that are taken in nuclear matter calculations always pertain to the 
sshell 4 ). The reason for adopting this correlation lies in its simplicity. At present the techniques 
for incorporating the realistic three-body AAA correlations for nuclear matter are not sufficiently 
developed as compared to those for the s-shell hypernuclei incorporated in this work. These affect 
the contribution from V^Jj N quite substantially, even to the extent of reversing its sign in the 
presence of the AAA correlations as can be seen in the present as well as in the \ 7 studies. The 
correct incorporation of the three-body correlations in nuclear matter may affect the results to 
quite an extent, particularly those at high densities [32]. The incorporation of these correlations is 
indeed a challenging task and is very much needed for the present work as well as for other related 
studies. 



§5. Conclusion and Comments 

From the results discussed in the previous section we note that the values of the spin-average 
strength V= 6.20 and 6.15 MeV give a reasonable description of the s-shell hypernuclei. We 
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have been able to provide a consistent account of the B\ values of \H, \H* , and %He using 
the realistic Argonne Vi$ NN interaction and Urbana IX NNN interaction alongwith A-nuclear 
interactions with appropriate correlations. Our results for B\ show very similar trends with the 
spin-average strength V of the AN interaction, and the ANN interaction parameters, C p and W a 
for all the s-shell hypernuclei considered. 

An important conclusion of our study is that about 25% to 50% of the + — 1 + splitting 
energy between the \H and \H * comes from the A A spin-dependent strength V a . The earlier 
studies [4,28,29] attribute a larger part of the splitting to the spin dependence of the two-body AN 
interactions. In contrast, our study indicates that the major part (~ 50% to 75%) of the splitting 
is generated by the three-body AAA forces. 

Our study on A-binding to nuclear matter shows that A is underbound. This indicates the 
fact that there is a need to include the three-body correlations while treating nuclear matter. 
This would require a different technique altogether and is a challenging problem in itself. Our 
analysis on j^O also indicate the importance of the non-central correlations. It is possible that 
the inconsistency between the results of Rcf. [16] and our s-shell results is due to their neglecting 
terms with higher than four-baryon clusters and thereby neglecting the contributions from the 
non-central clusters. Moreover, our values of the AAA interaction parameters, C p and W a are 
higher than those of Ref. [16] and which, in turn, would induce stronger AAA correlations. 

Contrary to the findings by Bando et al. [33] and Shinmura et al. [34] regarding the effect of 
tensor forces on the overbinding problem, we find that the tensor forces do not play a significant 
role. Further, separate studies by Carlson[28] and Hiyama et al. [29] and on four- and five-body 
hypernuclei have shown that the binding energies and the splitting energies are not reproduced 
correctly using the Nijmegen interactions which have strong tensor terms. However, the small 
suppression effects expected from AA tensor forces are already implicitly included in our phe- 
nomenological dispersive ANN force. Moreover, the Argonne v\g potential used in our study has 
a weak tensor part and this in fact, provides a much better binding to nuclei 5 '. 

ACKNOWLEDGMENTS 

We would like to thank Dr. R. B. Wiringa for providing his computer routines for the mass 3 
and mass 4 nuclei. We gratefully acknowledge Professor A. R. Bodmer for fruitful discussions and 
suggestions during the course of this work. 

REFERENCES 

[1] Dean Halderson, Phys. Rev. C61, 034001(2000) 
[2] Y. Yamamoto and H. Bando, Prog. Thco. Phys. Supp. 81, 9 (1985) 
[3] A. R. Bodmer, Q. N. Usmani, and J. Carlson, Phys. Rev. C29, 684 (1984) 
[4] A. R. Bodmer and Q. N. Usmani, Nucl. Phys. A477, 621 (1988) 

11 



[5] R. B. Wiringa, V. G. J. Stoks, and R. Schiavilla, Phys. Rev. C51, 38 
(1995) 

[6] J. Carlson, V. R. Pandharipande, and R. B. Wiringa, Nucl. Phys. A401, 
59 (1983) 

[7] B. S. Pudlincr, V. R. Pandharipande, J. Carlson, and R. B. Wiringa, 

Phys. Rev. Lett. 74, 4396 (1995) 
[8] Q. N. Usmani, Nucl. Phys. A340, 397 (1980) 

[9] R. B. Wiringa, R. A. Smith, and T. L. Ainsworth, Phys. Rev. C29, 1207 
(1984) 

[10] J. R. Bergervoet, P. C. van Campen, R. A. M. Klomp, J. L. de Kok, T. 
A. Rijkcn, V. G. J. Stoks, and J. J. de Swart, Phys. Rev. C41, 1435 
(1990) 

[11] V. G. J. Stoks, R. A. M. Klomp, M. C. M. Rentmeester, and J. J. de 

Swart, Phys. Rev. C48, 792 (1993) 
[12] J. Fujita and H. Miyazawa, Prog. Thcor. Phys. 17, 360 (1957) 
[13] I. E. Lagaris and V. R. Pandharipande, Nucl. Phys. A359, 331 (1981) 
[14] Thomas Schimert, D. J. Stubeda, M. Lemere, and Y. C. Tang, Nucl. 

Phys. A343, 429 (1980); Wang Xi-chang, H. Takaki, and H. Bando, 

Prog. Thcor. Phys. 76, 865 (1986) 
[15] R. K. Bhaduri, B. A. Loiseau, and Y. Nogami, Annals of Phys. 44, 57 

(1967) 

[16] A. A. Usmani, S. C. Piepcr, and Q. N. Usmani, Phys. Rev. C51, 2347 
(1995) 

[17] M. K. Pal, Theory of Nuclear Structure (Affiliated 

East- West, New Delhi) Appendix A, pg 606 
[18] R. B. Wiringa, Phys. Rev. C43, 1585 (1991) 

[19] A. Arriaga, V. R. Pandharipande, and R. B. Wiringa, Phys. Rev. C52, 
2362 (1995) 

[20] S. Murali, Ph.D. Thesis, Jamia Millia Islamia, New Delhi, India (1995) 
[21] J. Lomnitz-Adler, V. R. Pandharipande, and R. A. Smith, Nucl. Phys. 
A361, 399 (1981) 

[22] J. Carlson and R. B. Wiringa, in Computational Nuclear Physics, cd. by 
S. E. Koonin, K. Langankc, J. A. Maruhn, and M. R. Zirnbauer (Springer- 

12 



Verlag, Berlin, 1991) 
[23] N. Metropolis, A. W. Rosenbluth, M. N. Roscnbluth, A. H. Teller, and 

E. Teller, J. Chem. Phys. 21, 1087 (1953) 
[24] J. L. Forest, V. R. Pandharipande, and A. Arriaga, Phys. Rev. C60, 

014002 (1999) 

[25] Q. N. Usmani and A. R. Bodmer (private communication) 

[26] J. Pniewski and D. Zicminska, in Proceedings of the Conference on Kaon- 

nuclear interactions and Hypernuclei, Zvenigorod, pg 33 (1977) 
[27] M. Bedjidian, E. Descroix, J. Y. Grossiord, A. Guichard, M. Gusakow, 

M. Jacquin, M. J. Kudla, H. Piekarz, J. Pickarz, J. R. Pizzi, and J. 

Pniewski, Phys. Lett. 83B, 252 (1979), and references therein. 
[28] J. Carlson, in LAMPF Workshop on (n, K) Physics, ed. B. F. Gibson, 

W. R. Gibbs, and M. B. Johnson, AIP Conf. Proc. No.224 (AIP, New 

York, 1991) 

[29] E. Hiyama, M. Kamimura, T. Motoba, T. Yamada, and Y. Yamamoto, 

Nucl. Phys. A639, 169c (1998) 
[30] Q. N. Usmani and A. R. Bodmer, Nucl. Phys A639, 147c (1998) 
[31] Q. N. Usmani and A. R. Bodmer, Phys. Rev. C60 055215 (1999) 
[32] A. Akmal, V. R. Pandharipande, D. G. Ravenhall, Phys. Rev. C58 1804 

(1998) and references therein. 
[33] H. Bando and I. Shimodaya, Prog. Theor. Phys. 63, 1812 (1986) 
[34] S. Shinmura, Y. Akaishi, and H. Tanaka, Prog. Theor. Phys. 65, 1290 

(1981) 



13 



TABLE 1. AiV Correlations Parameters 





V 


KAN 


a\N 


Can 


Ran 






A 


6.20 


0.117 


0.50 


2.0 


1.0 


0.965 


0.80 




6.15 


0.110 


0.50 


2.0 


1.0 


0.970 


0.95 




6.10 


0.095 


0.50 


2.0 


1.0 


0.940 


0.95 


\H 


6.20 


0.12 


0.70 


2.0 


1.0 


0.95 


0.70 




6.15 


0.10 


0.70 


2.0 


1.0 


0.95 


1.20 




6.10 


0.08 


0.70 


2.0 


1.0 


0.95 


0.70 




6.20 


0.095 


0.70 


2.0 


1.0 


1.0 


0.70 




6.15 


0.065 


0.70 


2.0 


1.0 


1.0 


0.70 




6.10 


0.050 


0.70 


2.0 


1.0 


1.0 


0.70 



f Eq.(3.8); For all other correlation parameters refer to 16) 



TABLE 2. ANN Correlation Parameters for A He, %H, and \H* 



Parameters 


V=6.20 MeV 


y=6.15 MeV 


^=6.10 MeV 




0.364104 


0.311733 


0.257894 




0.006096 


0.004845 


0.003766 



f Eq.(3.11) 



TABLE 3. Variational results for 3 H and 4 He 



Components 


3 H 




Kinetic energy 


50.76(4) 


106.85(6) 


NN Potential energy 


-57.95(4) 


-129.30(6) 


NNN Potential energy 


-1.13(3) 


-5.27(7) 


Total energy 


-8.32(2) 


-27.71(6) 


r.m.s. (proton) 


1.585(3) 


1.478(2) 


r.m.s. (neutron) 


1.731(4) 


1.478(2) 


d-state probability 


0.0933(1) 


0.1512(2) 


All Energies are in MeV and Radii in fm 
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TABLE 4. Energy expectation values (calculated and fitted) 
for ^(Ground State) with V=6.20 MeV 



1/ 
V<J 


Wo 


c 


< i^cal > 


< tiftt > 


0.17 








10.58(05) 


10.60 


0.17 





1.0 


11.61(07) 


11.63 


0.17 


0.01 





9.97(05) 


9.97 


0.17 


0.01 


1.0 


10.96(05) 


10.90 


0.17 


0.02 





9.25(04) 


9.31 


0.17 


0.02 


1.0 


10.21(05) 


10.13 


0.17 


0.02 


2.0 


12.28(07) 


12.31 


0.17 


0.005 





10.23(05) 


10.29 


0.17 


0.005 


1.0 


11.25(05) 


11.27 


0.17 


0.005 


2.0 


13.58(08) 


13.60 


0.17 


0.015 


1.0 


10.50(05) 


10.52 


0.22 








10.82(05) 


10.71 


0.22 





1.0 


11.66(05) 


11.74 


0.22 


0.01 





10.02(05) 


10.09 


0.22 


0.01 


1.0 


10.90(05) 


11.01 


0.22 


0.01 


2.0 


13.36(08) 


13.30 


0.22 


0.02 





9.43(05) 


9.43 


0.22 


0.02 


1.0 


10.21(05) 


10.25 


0.22 


0.02 


2.0 


12.35(06) 


12.43 


0.22 


0.005 


1.0 


11.47(05) 


11.38 


0.22 


0.015 


1.0 


10.68(05) 


10.63 


0.27 








10.79(05) 


10.83 


0.27 





1.0 


11.79(06) 


11.86 


0.27 





2.0 


14.32(08) 


14.24 


0.27 


0.01 





10.26(05) 


10.20 


0.27 


0.01 


1.0 


11.13(05) 


11.12 


0.27 


0.01 


2.0 


13.41(08) 


13.41 


0.27 


0.02 





9.60(05) 


9.55 


0.27 


0.02 


1.0 


10.39(05) 


10.36 


0.27 


0.02 


2.0 


12.56(06) 


12.54 


0.27 


0.005 





10.50(05) 


10.52 


0.27 


0.005 


1.0 


11.46(05) 


11.49 


0.27 


0.015 


1.0 


10.81(05) 


10.75 



e =0.24, k A jv =0.12, a s =0.95, a a =0.7 
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TABLE 5. B\ as a function of coefficients y\_ 


-6 


iz V Vl 


V2 V3 


Vi 2/5 


ye 


\He 6.20 -0.92 


-285.51 1.28 


1701.78 1.50 - 


-37.62 


6.15 1.19 


-280.03 0.99 


1134.41 1.35 


-7.44 


6.10 -0.20 


-210.02 0.88 


2997.79 1.19 


A A HQ 

-44. Uo 


%H 6.20 2.27 


-61.23 0.35 


-141.55 0.68 - 


■10.53 


6.15 2.21 


-59.03 0.32 


453.47 0.52 


-8.14 


6.10 1.55 


-43.48 0.23 


298.23 0.37 


-4.96 


\H* 6.20 -0.89 


-112.62 0.37 


639.67 0.68 


-8.51 


6.15 -0.98 


-65.63 0.27 


208.58 0.45 


-6.68 


6.10 -0.81 


-53.21 0.12 


334.70 0.36 


-5^69 


t Includes contribution due to A-nuclear correlations 




TABLE 6(A). 


Variational Results for \He 




Components 


V = 6.20 


V = 6.15 


V = 6.10 


Nuclear Kinetic energy^ 


128.38(74) 


128.16(71) 


125.35(69) 


NN Potential energy^ 


-139.97(73) 


-140.27(69) 


-139.31(67) 


NNN Potential energy 


-5.97(8) 


-5.91(8) 


-5.73(8) 


A Kinetic energy 


11.64(15) 


11.11(15) 


9.13(13) 


AN P.E (central) 


-23.65(32) 


-21.69(30) 


-17.55(27) 


AJVP.E (spin) 


-0.0138(1) 


-0.0311(3) 


-0.0281(3) 


AN Space exch. contribution 


0.763(13) 


0.578(10) 


0.395(7) 


ANN P.E (Total) 


-1.91(9) 


-2.67(10) 


-3.53(10) 


ANN P.E (T.P.E.) 


-8.43(16) 


-8.77(17) 


-8.20(17) 


ANN P.E (dispersive) 


6.52(11) 


6.10(11) 


4.67(10) 


Total Energy 


-30.75(14) 


-30.76(18) 


-31.27(12) 


B A 


3.03(15) 


3.05(19) 


3.56(13) 


r.m.s. radius (proton) 


1.376(2) 


1.379(2) 


1.389(2) 


r.m.s. radius (neutron) 


1.377(2) 


1.379(2) 


1.389(2) 


d-state probability 


0.1568(2) 


0.1568(2) 


0.1557(2) 



f Includes contribution due to A-nuclear correlations 
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TABLE 6(B). Variational results for \H 


Components 


V = 6.20 


V = 6.15 


V = 6.10 


Nuclear Kinetic energy' 


65.76(47) 


62.82(47) 


60.34(46) 


A//\/ T-Vitr»nti^l niKifffui 
IV IV JT ULt^llLldjl H^y 


fifi Q1 (&7\ 




UO .UUl 4U 1 


A//\//\/ T-^f~\t HTlt 1 'A 1 OTlPriTU 
IV IV IV JT U LL.I1L1C11 L.llt^lJiV 


l.OOl O 1 


1 95^^ 

J. .ZjO 1 O J 




A Kinetic energy 


7.17(11) 


5.99(9) 


4.77(9) 


AN P.E (central) 


-14.00(22) 


-11.55(20) 


-8.84(18) 


AA P.E (spin) 


-0.291(3) 


-0.363(4) 


-0.300(4) 


AN Space exch. contribution 


0.333(9) 


0.221(6) 


0.143(4) 


AAA P.E (Total) 


-1.33(6) 


-1.36(6) 


-1.48(6) 


A AW P.E (T.P.E.) 


-2.87(8) 


-2.63(9) 


-2.47(9) 


ANN P.E (dispersive) 


1.54(4) 


1.27(4) 


0.99(3) 


Total Energy 


-10.61(6) 


-10.59(5) 


-10.16(5) 


B A 


2.28(6) 


2.27(6) 


1.83(6) 


r.m.s. radius (proton) 


1.408(2) 


1.435(3) 


1.460(3) 


r.m.s. radius (neutron) 


1.516(3) 


1.547(3) 


1.577(3) 


d-state probability 


0.0984(1) 


0.0960(1) 


0.0962(1) 



f Includes contribution due to A-nuclear correlations 



TABLE 6(C). 


. Variational results for \H* 




Components 


V = 6.20 


V = 6.15 


V = 6.10 


Nuclear Kinetic energy' 


63.78(49) 


60.58(47) 


57.96(46) 


NN Potential energy^ 


-65.71(48) 


-63.64(46) 


-62.36(45) 


NNN Potential energy 


-1.34(3) 


-1.27(3) 


-1.20(3) 


A Kinetic energy 


6.40(12) 


4.78(9) 


3.4516(8) 


AA P.E (central) 


-12.35(25) 


-9.42(19) 


-6.56(17) 


AAA P.E (spin) 


0.086(1) 


0.083(1) 


0.065(1) 


AA Space exch. contribution 


0.311(8) 


0.202(6) 


0.110(4) 


AAA P.E (Total) 


-0.57(6) 


-0.62(5) 


-0.73(5) 


AAA P.E (T.P.E.) 


-3.11(9) 


-2.32(8) 


-1.98(8) 


AAA P.E (dispersive) 


2.54(7) 


1.70(5) 


1.24(5) 


Total Energy 


-9.40(7) 


-9.31(5) 


-9.26(5) 


B A 


1.08(7) 


0.99(5) 


0.94(5) 


r.m.s. radius (proton) 


1.431(3) 


1.467(3) 


1.494(3) 


r.m.s. radius (neutron) 


1.542(3) 


1.585(3) 


1.618(3) 


d-state probability 


0.0995(1) 


0.0980(1) 


0.0969(1) 



f Includes contribution due to A-nuclear correlations 
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TABLE 7. Breakup of the + — 1 + splitting contributions 



Contribution 


V = 6.20 


V = 6.15 


V = 6.10 


v a 


0.377(3) 


0.446(4) 


0.365(4) 


Vann 


0.76(8) 


0.74(8) 


0.75(8) 


Total 


1.137(80) 


1.186(80) 


1.115(89) 


Energy differences 


1.21(9) 


1.28(7) 


0.90(7) 



The first row gives the contribution to splitting from V CT . The second row gives the 
contribution arising from Va AT N • The third row gives the total of V CT and VAN N ■ 
The last row gives the actual calculated energy difference between \H and \H* . 



TABLE 8. Results for Nuclear Matter calculations 



V 


Po 


-D 


<T A + V AN > 


Spc. exch. 


Vann 


6.20 


0.162 


-21.525 


-77.617 


8.561 


47.530 


6.15 


0.162 


-17.819 


-72.628 


6.644 


48.164 


6.10 


0.162 


-11.727 


-67.729 


4.798 


51.205 



All energies are in MeV. p Q is the normal nuclear matter density in fm 
The third column gives the well depth D. The fourth column gives the value 
of D without the ANN forces and the space-exchange contribution, i.e. for 
6 = 0. The fifth column gives the reduction in the contribution to D due to 
the space-exchange part (Spc. exch.) of the AN potential. The last column 
gives the contribution due to the three-body ANN forces. 
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Figp 2 . B ^ as a ia nclioti of W 




